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Gradient expansion approach to nonlinear superhorizon perturbations II 

— a single scalar field - 

Yoshiharu TANAKA* and Misao SASAKP 
Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

We formulate nonlinear perturbations of a scalar field dominated universe on super-horizon scales. 
We consider the case of a single scalar field. We take the gradient expansion approach. We adopt 
the uniform Hubble slicing and derive the general solution valid to 0(e 2 ), where e is the expansion 
parameter associated with a spatial derivative, which includes both the scalar and tensor modes. 
In particular, the 0(e 2 ) correction terms to the nonlinear curvature perturbation, which become 
important in models with a non-slowroll stage during inflation, are explicitly obtained. 
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I. INTRODUCTION 



The cosmic microwave background (CMB) anisotropics recently observed by WMAP strongly support the infla- 
tionary cosmology, and theoretical predictions of various models of inflation seem to be well estimated by linear 
perturbation theory, with the cosmological perturbations generated from quantum vacuum fluctuations which are 
well approximated by Gaussian random fields Nevertheless, there has been a growing interest in detection of 
possible derivations from the Gaussian statistics. It was suggested that a deviation from the Gaussian stastics may be 
used to distinguish models of inflation Q, and it may indeed be detected by PLANCK 0] in the near future. Thus, 
making clear predictions on the non-Gaussianity from inflation have become one of the urgent issues of the inflation- 
ary cosmology. Since the nonlinearity is essential for the generation of non-Gaussian perturbations, it is necessary to 
develop a nonlinear cosmological perturbation theory to evaluate the non-Gaussianity from inflation. 
, Our goal is to formulate a nonlinear theory with which we can calculate non-Gaussianities from any models of 
■ inflation. The traditional approach is to develop a second-order perturbation theory [J, |5|, [6|, [8|, [9( . However, here we 
OO | adopt a different one, the gradient expansion approach, in which nonlinear perturbations are solved by invoking spatial 
gradient expansion under the assumption that spatial derivatives are small compared to time derivatives. Technically, 
we introduce an expansion parameter, e, and associate it with each spatial derivative. Then we expand field equations 
in terms of e and solve them order by order in e iteratively. The gradient expansion a ppr oach has been developed and 
studied by many authors previously [O, [H [H Q M, M, _3_ M, El S [ll 111 111, 13, [11, [11, [11 . 

In cosmological situations, the gradient expansion approach is valid on scales greater than the Hubble horizon scales, 
and an advantage of the approach is that we can calculate perturbations to full nonlinear order in their amplitudes. 
This is particularly useful when dealing with the general non-Gaussianity for which it may be necessary to evaluate 
not only second-order perturbative corrections but also higher order corrections. 

At the leading order in the gradient expansion, i.e. neglecting all the spatial gradients, Lyth, Malik and Sasaki fi~fl ] 
studied nonlinear scalar curvature perturbations, proved the nonlinear SN formula, and constructed a gauge invari- 
ant (time-slice independent) nonlinear scalar curvature perturbation. Although this leading order approximation is 
sufficient for a large class of inflation models, there are models for which it is necessary to take into account the next 
order corrections. For example, in the context of the standard linear perturbation theory, Leach et al. pointed out 
that there can be enhancement of the comoving curvature perturbation on superhorizon scales, where it is usually 
conserved, even in a single field model if the slow-roll condition is temporarily violated [_8l |. There it was shown that 
0(fc 2 ) corrections to the curvature perturbation on superhorizon, where k is the comoving wavenumber, play a crucial 
role for the enhancement. Because 0(k 2 ) corrections correspond to 0(e 2 ) terms in gradient expansion, this implies 
that it is necessary to include 0(e 2 ) terms in such a model. Then, we expect that the enhancement may give rise to 
large non-Gaussianity. Indeed, Chen et al. [lfl ] numerically found in a single-field inflation that large non-Gaussianity 
can be generated if the slow-roll condition is temporarily violated, using third-order action derived by Maldacena [![. 

In this paper, focusing on a single-field inflation, we formulate the gradient expansion to 0(e 2 ) on the uniform Hubble 
slicing. In most of the previous studies, either only the leading order terms in gradient expansion was discussed or 
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the choice of time-slicing was not quite adequate for the study of non-Gaussianity from inflation. Here we adopt the 
uniform Hubble slicing because the curvature perturbation on this slicing directly determines the initial condition for 
the CMB anisotropies and the large scale structure of the universe. 

We employ the (3 + l)-decomposition of the Einstein equations, and consider a single scalar field with an arbit rary 
potential. We then derive the general solution for all the variables. As discussed in the case of a perfect fluid in [31j , 
we find that the identification of the tensor mode in the spatial metric is rather arbitrary, depending on how one fixes 
the spatial coordinates, as a reflection of general covariance, while it can be unambiguously identified in the extrinsic 
curvature of the metric. 

This paper is organized as follows. In Section [H] we define basic variables, and describe the assumptions we adopt 
in gradient expansion. Then we present the general solution for all the physical quantities to 0(e 2 ) on the uniform 
Hubble slicing. In Section Hill we discuss the validity of the assumption adopted in Section [II] by appealing to linear 
theory and by considering the vacuum fluctuations of the scalar field at and around horizon crossing. We conclude the 
paper in Section IV. In Appendix [Al the basic equations are presented. In Appendix [B] the estimation of the orders 
of physical quantities in powers of e is given. In Appendix \C\ the general solutions for all the physical quantities are 
derived. 



II. GRADIENT EXPANSION 



A. Basic variables 



In the (3 + l)-decomposition, the metric is expressed as 

ds 2 — gfurdx^dx" 

= {-a 2 + 8 k 8 k )dt 2 +20 i dx i dt + -f ij dx i dx j , (2.1) 

where a, B % (3 l = 7 y /3j), and 7^ are the lapse function, shift vector, and the spatial metric, respectively. We rewrite 

lij as 

lij {t,x k ) = a 2 {t)^{t,x k )%{t,x k ); det(7y) = l, (2.2) 

where the function a(t) is the scale factor of a fiducial homogeneous and isotropic background universe. 
The extrinsic curvature is defined by 

Kij = -Vii^ , (2.3) 
where n M = (—a, 0, 0, 0) is the vector unit normal to the time slices. We decompose the extrinsic curvature as 



Kij = ^fK + ip A a 2 A l0 ; K = K %3 , (2.4) 



where Aij represents the traceless part of Kij . The indices of K^ are to be raised or lowered by 7 U and 7^ , and the 
indices of by 7 y and 7^ . 

The stress-energy tensor for a single scalar field is 

V = V^V^ - ~<v(V a 0V a + 21/(0)). (2.5) 

We define the local Hubble parameter as 1/3 of the expansion of the unit normal vector n M , which is equal to —1/3 
of the trace of the extrinsic curvature in our convention, 

3H = -K = h 6^- — , (2.6) 

aa oop a 

where a dot ' denotes d/dt and Di is the covariant derivative with respect to 7ij. In the following, we adopt the 
uniform Hubble slicing. For this slicing, we have 



and Eq. (|2.6p implies 



lf(t) = ~, (2.7) 



Q = 1 + W^w- (2 - 8) 
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B. Expansion scheme 



We investigate nonlinear superhorizon perturbations with the gradient expansion approach, which is called by 
various names by various authors, the quasi- isotropic expansion fl2l. [T3, [l9l. u)\ . the anti-Newtonian approximation 
0, , the spatial gradient expansion 0, El, H3, Hlf, [H, [H, [25| , or the long wavelength approximation [26|, [2i| . 
In this approach, we assume that the characteristic length scale L of inhomogeneities is always much larger than 
the Hubble horizon scale, L 3> H^ 1 ~ t. We introduce a small parameter e, and assume that L is of 0(l/e). 
This assumption is equivalent to assuming that the magnitude of spatial gradients is given by ditp = ip x O(e), 
did — ax O(e), etc.. In the limit L — » oo, i.e., e — * 0, the universe looks locally like a FLRW spacetime, where 
'locally' means as seen on the scale of the Hubble horizon volume. It is noted that physical quantities which are 
approximately homogeneous on each Hubble horizon scale can vary nonlinearly on very large scales. 

The local homogeneity and isotropy imply that (i % — 0(e) and dt'jij = 0(e) because the local FLRW equations 
should be realized in the limit e — > 0. However, we further assume that (3 Z = 0(e 3 ) and 9t7y = 0(e 2 ). Technically, 
these additional assumptions make the analysis of 0(e 2 ) correction terms much simpler, as discussed in Appendix [Cll. 
Physically, of course, it is necessary to justify them. For the former assumption on (3 l , since it is just a matter of 
choice of the spatial coordinates, this does not cause any loss of generality. In fact, once we obtain the solution, 
it is straightforward to express it in a more general spatial coordinate system by the coordinate transformation 
x l — ► x l — F l (t,x k ) such that dtF 1 = 0(e), corresponding to /3 l = 0(e) in the new coordinate system. As for the 
latter assumption on dt'jij, however, it is not simply a matter of choice. In Sec. IIII| with the help of the result from 
linear theory, we give a convincing argument, if not rigorous, that this assumption is indeed satisfied for perturbations 
arising from the vacuum fluctuations. To summarize, our basic assumptions are 

/3 4 = 0(e 3 ) , daij = 0(e 2 ) . (2.9) 

Applying these assumptions to the Einstein-scalar field equations, we find 

V; = 0(1), a-l = 0(e 2 ), c^ = 0(e 2 ), A tJ = 0(e 2 ) . (2.10) 

These estimates are derived in Appendix [B"l Here one comment is in order. The fact that dtip = Q(e 2 ) means ip is 
conserved if the 0(e 2 ) corrections can be neglected. This result was derived by Salopek and Bond [1J| for a single 
scalar field system, and by Lyth, Malik, and Sasaki [ll[ for more general systems. 



C. General solution 



Here we present the general solution for all the physical quantities, valid to 0(e 2 ) in gradient expansion. We defer 
the derivation to Appendix [U] because it is not much different from the one we gave in the previous paper [3l[, except 
for the fact that the present paper deals with the case of a scalar field while the previous paper dealt with a perfect 
fluid. 

The general solution is 



= 1 



a 3 (f> 



in 

An 



{2) C{x k ) (^+^ / a(t')dt') + ( 2 ) D{x k )^ I 
(0) L(x k )(l + - J (a-l)Hdt') 
m fik(x i )(5$-2 (2) F k j (x i ) 



a(t')dt' 



dt' 

(2)C l] {x k ) 



a(t")dt" -2 {2) C k j{x e ) 



dt' 

. ^0 



(j) (t) + {2) C(x k )^ 1 / a (t')dt' + (2) D{x k )^- 



(2.11) 
(2.12) 
(2.13) 
(2.14) 
(2.15) 



where the index (n) is attached to a quantity of 0(e") except for the scalar field. For notational simplicity, the full 
scalar field is denoted by 0, and the lowest order scalar field im<j>, which depends only on time, is denoted simply by 
4>. The time £» is an arbitrary reference time and </>* = cj>(t*). The tensor (a)-Pij and the scalar r 2 \C are given by 
Eqs. (|C14|) and (|C34[) , respectively, as functions of ( )£ and (o)fij- The te nsor ( 2 )Cjj is traceless with respect to (o)/i?, 
The function is related to ( 2 )Cij through the momentum constraint (|C36[) . 
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To clarify the physical role of these freely specifiable functions, let us first count the degrees of freedom. Since (2)C 
and (2)Fij are determined by and (o)/y, they have no degree of freedom. Since the determinant of (o)/y is unity, 
it has 5 degrees of freedom, and since (2)CV, is traceless, it also has 5 degrees of freedom. In addition we have ( )L and 
(2)D, each of which counts 1 degree of freedom. The momentum constraints which consist of 3 equations, relate (2)CV, 
to (2)D, and reduce the total degrees of freedom by 3. So, adding together, the total number is 5 + 5 + 1 + 1 — 3 = 9, 
while the true physical degrees of freedom are 2 + 2 + 2 = 6, where 2 are of the single scalar field (1 for a growing 
mode and 1 for a decaying mode) and 2 + 2 are of the gravitational waves (2 for the metric and 2 for the extrinsic 
curvature). This implies that there still remains 3 degrees of freedom. 

As discussed in the case of a single perfect fluid in the previous paper [3l|, the remaining 3 degrees of freedom 
come from the spatial covariance, that is, from the gauge freedom of purely spatial coordinate transformations 
x l — > x l = f l (x j ). Thus we may regard that (o)/ij contains these 3 gauge degrees of freedom. The correspondence 
between the nonlinear solutions and the solutions in linear theory can be understood from the time dependence of the 
nonlinear solutions. So, we see that ( )£ an( i (o)fij represent growing modes, and (2)-D and (2)Cjj represent decaying 
modes. 

To summarize, the degrees of freedom contained in the freely specifiable functions can be interpreted as 

( )i • • • 1 = 1 (scalar growing mode) , 

(o) fij ' ' ' 5 = 2 (tensor growing modes) + 3 (gauge modes) , 

(2) Cij ■ ■ ■ 5 = 2 (tensor decaying modes) + 3 (constraints) , 

( 2 )D ■■■ 1 = 1 (scalar decaying mode) . (2-16) 

To fix the gauge completely, one has to impose 3 spatial gauge conditions on (o)/y to extract the physical tensor 
degrees of freedom. As discussed in Ref. 31], this cannot be done in a spatially covariant way because (o)/ij is the 
metric and any covariant derivative of it vanishes identically. Thus the tensor modes cannot be extracted out from 
(o)/ij unless we introduce a certain 'background' metric. This is an important difference from the linear case in which 
there exists a background metric. 

In contrast, as for the extrinsic curvature, we may identify the tensor modes in it, because its transverse-traceless 
part can be extracted unambiguously [30l l3ll|. Namely, the transverse part of ^L 6 Aij can be determined uniquely 
and it can be identified as the tensor modes. As noted in Ref. [3lj |. this implies that the tensor modes in the extrinsic 
curvature are determined non-locally, and they exist even for a trivial ( ) fij , say for ( ) f t j = Sij . This generation of 
tensor modes in the extrinsic curvature is a result of nonlinear interactions of the scalar modes. We note, however, 
that it is not obvious if the tensor modes we identified can be called gravitational waves. To make a clear connection 
between the tensor modes and gravitational waves, it is necessary to evolve the system until the scale of interest 
becomes sufficiently smaller than the Hubble scale. But this is beyond the scope of the present paper. 

III. VALIDITY OF THE ASSUMPTION IN THE LINEAR LIMIT 

In this section, we discuss the validity of our central assumption, dt"fij = 0(e 2 ). With the help of linear theory, we 
argue that it can be physically justified. We first consider the equation for the curvature perturbation during inflation, 
and argue that the condition dtjij = 0(e 2 ) is naturally satisfied. Then considering the quantum fluctuations as the 
source of the curvature perturbation, we explicitly show that this assumption holds not only in the case of slow-roll 
inflation but also in the case of the Starobinsky model, in which the slow-roll condition is temporarily violated due 
to a sudden change in the slope of the potential. 

Since gradient expansion is effective only on superhorizon scales, as the initial condition for the quantities we 
calculate, we need to know their behavior when their wavelength crosses the Hubble horizon radius. To do so, we 
assume that the linear perturbation is a sufficiently good approximation up to and around the Hubble horizon scale. 
For linearized quantities, we follow the notation of Kodama and Sasaki [2]. 

As usual, we expand the linearized quantities in terms of spatial harmonics. The background is a spatially flat 
FLRW universe. We introduce the scalar, vector and tensor harmonics, Y^ 1 ) and Y^)-, respectively, 

(A + k 2 )Y k = , 

( A + fc 2 )y«= ; 8^=0, 

(A + k 2 )Y^ = ; 5*Y hij = = , (3.1) 

where A = d l di is the 3-dimensional Laplacian. For the sake of notational simplicity, we suppress the mode index 
k below. For a universe dominated by a scalar field, it is known that the vector perturbations are not be excited. 
Therefore, we ignore the vector modes. 
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We express the metric as 

ds 2 = a{r,) 2 [-(1 + 2AY)dif - 2BY j dx j dr] + (<% + 2H L Y8 lJ + 2H T Y lj + 2H^ ) Y l f ) )dx t dx j , (3.2) 
where dr\ = dt/a{t) is the conformal time, and 



-2 



didj - ^.\,A 



Y. (3.3) 



The correspondences of these to the metric components defined in Section |TT] are 

a = l + AY, pj^-aBYj, ij 2 = l + H L Y, % = % + 2H T Y ij + 2H^ ) Y§ ) . (3.4) 

In passing, we note that our assumption 9t7y = 0(e 2 ) corresponds to dtHx = 0(e 2 ) and dtH^ = 0(e 2 ) in the linear 
limit. We also introduce the quantities 

K = H L + ^, (3.5) 
a g = k- r H' T -B, (3.6) 
K = -A + Afi + H- X H' L = -A + H-'W - A(7 S , (3.7) 

where a prime ( ') denotes a conformal time derivative djdr\. These quantities are known to be independent of the 
choice of the spatial coordinates but depend only on the choice of time-slicing. The first one, 1Z is called the curvature 
perturbation because the spatial curvature scalar R[j] is given by 

i?[ 7 ] = ~4 A ^ y ] • ( 3 - 8 ) 
The variables K, and a g represent the perturbations in the extrinsic curvature, 

K = —3H (1 + ICY) , A i:j = --cTgYij . (3.9) 

Under the gauge transformation induced by an infinitesimal change of time-slicing fj — r\ + TY, 1Z, a g and K, 
transform as 

K = K-HT, (3.10) 

°g = cr 9 - fcT, (3-11) 

K, = 1C + H- 1 (n 2 -W + ^t)t , (3.12) 



where TL = aH is the conformal Hubble parameter. 

Since the standard linear calculation gives the curvature perturbation on the comoving slicing, we need to relate 
the quantities on the comoving slicing to those on the uniform Hubble slicing. In particular, analytic expressions for 
the comoving curvature perturbations in the Starobinsky model resulting from the quantum fluctuations are given 
in [27| . Here and in what follows, we denote a quantity on the comoving and uniform Hubble slicing by the subscripts 
c and H, respectively. Thus, we first express the geometrical quantities on the comoving slicing in terms of 1Z C . Then 
we consider an infinitesimal transformation from the comoving slicing to the uniform Hubble slicing. 

The (0, /i)-components of the Einstein equations on the comoving slicing give 

6G° = n 2 5T°; 

2[m 2 A c + Hk{a g ) c - mn' c - k 2 K c ] = n 2 (j)' 2 A c , (3.13) 
*G§ = k 2 5T?; 

HA c -TZ' c = 0, (3.14) 
where k 2 = 8irG. From these, we have 

A c = hn' c , k(a g ) c = 1 (k 2 K c + > ( 3 - 15 ) 
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where it is noted that the background equation gives the relation, 



a^ = _o(-V = W a -W / = ^' 2 . (3.16) 



Then the second equality of Eq. (|3.T[) gives 



Since K, is zero on the uniform Hubble slicing by definition, we have 

k 2 



Q = ICh =IC C +H- 1 (H 2 -H' + ?-)T. (3.18) 



This gives T in terms of /C c , which in turn is given by Eq. (I3.17j) . The result is 



Thus, we obtain 

3 ; H 3^0/2 + fc2 y 

(ka g ) H = [H' T - kB) H = - ^ 2 J + g (^fc 2 ^ c + -^-tt' c J . (3.21) 

Since -B = 0(A;) because f3 l — 0(e), the above results imply that our assumption dt^fij = 0(e 2 ), which corresponds to 
H' T = 0(k 2 )HT in the linear limit, is justified if 1Z' C = 0(k 2 )TZ c on superhorizon scales. 

Now let us argue that this is indeed so in general. We know that 1Z C satisfies the equation, 

„ z , o ad>' ad> 

K + 2-K + k 2 U c = ; z =-^ = -£. (3.22) 

We consider an inflationary universe in which we approximately have 7i ~ —1/V- There are two independent solutions 
u and v, which are functions of kr\. We may fix these solutions such that they behave in the superhorizon limit 
k/H ~ k\?]\ -> as 



= [l + 0(k 2 )] u ; uo — const. , 

v = [1 + 0(k 2 )] v ; v = % , (3.23) 

drj'/z 2 ^') 



Vk 

where rjk ~ —A: -1 is the horizon crossing time, Ti.(rjk) = k, and we have used the fact that the long wavelength 
expansion of Eq. (|3.22|) would give corrections only in powers of k 2 . Then the general solution for 1Z C is given by 

TZ C = Au + Bv , (3.24) 

where A and B are constants of order unity in the sense of the k expansion. Thus at fc|?7| <C 1, we have 

K c = Au a , n' c = O(k 2 )Au + Bv' Q . (3.25) 

Now, assuming that the violation of the slow-roll condition is not too strong, which is necessary for the spectrum to 
be approximately scale-invariant, we have z oc a oc ry^ 1 , hence 
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Thus we conclude that 1Z' C = 0(k 2 )lZ c on superhorizon scales. 

To reinforce the above argument, let us consider the Starobinsky model in which the slow-roll condition can be 
significantly violated [13] • The Starobinsky model has a sudden change in its slope at <j> = cpo such that 



V{<j>) = 



V + A+{cf) - <j) ) for <j)><j)o, 
Vo + A-(<t> - <t> ) for cj)<(j)o. 



(3.27) 



where A + , A- and <po are assumed to be positive so that the scalar field evolves from a large positive value of 
toward = 0. Then the background scalar field 4> satisfies 



+ {A+-A-)e 



-3H(t-t ) 



for 
for 



> 00 ■ 
< 00 ■ 



(3.28) 



where to is the time at which (j> — 4>o, and the de Sitter approximation, 3H 2 — k 2 Vq is assumed to be valid. Thus 
the scalar field slow-rolls at > 0o, and violates the slow-roll condition temporarily at < 0q- The evolution is 
decelerated if A + / A- > 1, or accelerated if A+/A- < 1, compared to the slow- roll evolution. 

The curvature perturbation on the comoving slicing during the non-slow-roll regime at < 0o is [13] 



iH 2 



V2k 3 / 2 
a(k) 



(a(fc)e- lfc "(l + ikrj) - P(k)e lkv (l - ik v )^j , 
k 



m = - 3 i ( — 



A_ 

1. 



k 

likjka ^0 



7.2 

1 + 



ikp 
k 



(3.29) 
(3.30) 
(3.31) 



where ko = (a(io)-ffo) -1 . On superhorizon scales, we have k\rj\ <C 1. The standard slow-roll case corresponds to 
A + = A-. In this case, we immediately see that a(k) — 1, f3(k) — 0, and 1Z' C — 0(k 2 )H c on superhorizon scales. 

We now turn to the case A + ^ A-. We assume A^/A + = 0(1) so that the slow-roll condition is not severely 
violated. For k < ko, we expand the three exponentials in Eq. (|3.29[) to obtain 



3iH 3 



V2k 3 / 2 A^ 



1 



A+ (k v f 



A 



1 + — - 1 1 + 2ko V 



[k r)) s 



5(M)' 



•0(A) 



Thus we have TZ' C = 0(k 2 )TZ c . For k > k , we have a(k) = 1 + O(k /k) and j3{k) = O(k /k). Hence 

2 



0(k 2 



O ((kr)) 2 



0(k 2 ) 



(3.32) 



(3.33) 



Since k > ko, we see the right-hand side of this is also of 0(k 2 )lZ c . Thus we have shown that 1Z' C = 0(k 2 )lZ c holds 
even for the Starobinsky model in which the slow-roll condition can be violated. This in turn implies H' T = 0(k 2 )HT- 



(2) r -. 

Next, we consider the tensor perturbation Hrp . The evolution equation is 



H. 



(2)// 



a 



k 2 Hi 2) =0. 



(3.34) 



This equation has the same structure as Eq. (|3.22[) for 1Z C if we replace z with a. Then we can repeat exactly the 

same argument and conclude that H^' — 0(k 2 )Hj? . 

Again, this conclusion can be supported by considering the quantum fluctuations explicitly. When the de Sitter 
approximation H =const. holds, the tensor perturbation from the quantum fluctuations during inflation are given by 



rr(2) _ 



2irH 
W 2 



(kr) 



-ikrj 



Taking the conformal time derivative of this equation, we obtain at k\r)\ <C 1, 



H 



(2)/ 



2irH 



ik 2 rf)e 



-ikrj 



(2) 



0{k 2 )H^ 



(3.35) 



(3.36) 



Thus, provided that the linear theory is a good approximation up to scales corresponding to a few e- foldings after 
horizon crossing, our assumption dtjij = 0(e 2 ) is justfied for perturbations originated from the quantum fluctuations 
inside the horizon. 
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IV. CONCLUSION 



In this paper, taking the gradient expansion approach, we have investigated nonlinear perturbations on superhorizon 
scales in a universe dominated by a single scalar field. We have derived the general solution for all the physical 
quantities valid to second order in the spatial gradients on the uniform Hubble slicing. In particular, an expression 
for the nonlinear curvature perturbation, which plays a central role in the evaluation of a possible non-Gaussianity 
from inflation, has been obtained. 

Parallel to our previous paper (3lj |. we have identified the tensor modes in the extrinsic curvature, while the 
identification of the tensor modes in the metric is arbitrary unless we specify the background spatial metric. This is 
an important difference from the case of linear theory in which the background metric is uniquely given. 

In our analysis, we have adopted a non-trivial assumption dtjij — 0(e 2 ), which cannot be justified within the 
context of gradient expansion. To justify it, we have appealed to linear theory, and argued that the assumption 
dtlij = 0(e 2 ) is naturally satisfied in the linear limit. In particular, as an explicit example, we have considered the 
Starobinsky model, in which the slow-roll condition is temporarily violated due to a sudden change in the slope of the 
potential, and we have explicitly shown that the quantum fluctuations indeed satisfy the assumption d t jij — 0(e 2 ). 

As mentioned in Introduction, in models in which the slow-roll condition is temporarily violated as in the case of 
the Starobinsky model, the 0(k 2 ) corrections to the curvature perturbation, which correspond to 0(e 2 ) corrections in 
gradient expansion, may play a crucial role in the determination of the final amplitude of the curvature perturbation. 
Since the result of this paper is valid to 0(e 2 ), it provides a very useful tool to investigate the nonlinear behavior 
of the curvature perturbation on superhorizon scales for a wide class of models including those which violate the 
slow-roll condition. In particular, matching our result with the quantum fluctuations on scales slightly beyond the 
horizon scale, the non-Gaussianity arising from the nonlinear dynamics of the scalar field on superhorizon scales can 
be studied. Applications to specific models of inflation will be considered in a future publication. 
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We employ the (3 + l)-formalism of the Einstein equations and then the dynamical variables are 7^ and K^. The 
(n, n) and (n,j) components of the Einstein equations give the Hamiltonian and momentum constraint equations, 
respectively, while the (i, j) components gives the evolution equations for Kij. The evolution equations for 7^ are 
given by the definitions of the extrinsic curvature (|2.3p . 

In the present case, the Hamiltonian and momentum constraints are 
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APPENDIX A: BASIC EQUATIONS 



Here we show the basic equations for nonlinear quantities. The Klein-Gordon equation is 




(Al) 




8ttGJ 7 - 



= 16itGE , 



(A2) 



(A3) 



E = T^nPrt 



T^ Yj ■ 



(A4) 



The evolution equations for 7^ are given as 



{d t - (3 k d k )^ + — xl> = ^{-aK + d k p k }, 
la 

{d t - p k d k )%j = -2aAij + j ik djf3 k + % k dil3 k - ^d k p k , 



(A5) 



(A6) 



where ' = d/ dt. 
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The evolution equations for K%j are given as 



(fit ~ l3 k d k )K = tt (AyA« + -K 2 ) - D k D k a + ^Ga{E + S k k ), (A7) 



-a(KA l3 - 2A lk A k j ) 



+A ik d j (3 k + A jk d^ k - ^A t3 d k p k - - 1 fS k k ), (A8) 

where Rij is the Ricci tensor of the metric 7^ , R = 7 y i2y , Dj is the covariant derivative with respect to 7^ , and 

Sy = T y , 5 fe fe = 7 fci 5 ifc . (A9) 

APPENDIX B: ORDER ESTIMATION 

Here we evaluate the order of magnitude of the basic variables using the equations presented in Appendix [A] and 
the assumptions given by Eq. (|2.9[) , namely 

(3* = 0(e 3 ) , 9t7y - 0(e 2 ) • (Bl) 

With these assumptions, Eqs. (|A3[) and (|A6[) on the uniform Hubble slicing yield 

J 3 = Q(e 3 ). (B2) 



Using Eq. (|2.5[) . J, is express as 



Jj = --d 4>d^ + —drfdjcf) . (B3) 



We expand <p as (j> — m)0 + (i)0 + (2)4> + ' ' ' • Then, to satisfy Eq. (|B2[) . (o)0 should depend only on time, and m0 
should vanish. Thus, we obtain 

(0)0 = (o)^(*) . (i)0 = O. (B4) 
For notational simplicity, we denote the full scalar field by <j) and (o)0 by 4> in the following. Thus 

4>(t,x l ) = 4>{t) + {2) cj>{t,x i ) + --- . (B5) 
From the O(e ) part of the Hamiltonian constraint, Eq. (|A2[> . we have 

iA- 2 (;)=3tf 2 = 8^G(-^ + U(0)), (B6) 

where we have also expanded a as (o) a + (2) a + ' ' ' ■ From this equation, we find that (q)Q depend only on time, 

( )Q = (0) a(t). (B7) 



Since is spatial homogeneous, we may choose the time coordinate to set it to unity, ( ) a = 1- Thus, from Eqs. (|A5|) 
and 3H = — K we have 

= -3- (2) a + ^(l- (2) a) + O(e 4 ). (B8) 

We see dtip — 0(e 2 ) from this equation. 

To summarize, the orders of magnitude of the basic metric quantities are 

V; = 0(1), /3 l = 0(e), a-l = 0(e 2 ), dt% = 0(e 2 ) , dfijj = 0(e 2 ) , l y - = 0(e 2 ) . (B9) 

Actually, as for when we solve the equations to 0(e 2 ) in Appendix [Cl we assume ft 1 = C(e 3 ). But since the choice 
of /3 l does not affect the temporal behavior, it does not affect the generality of the solution. 
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APPENDIX C: DERIVATION OF THE GENERAL SOLUTION 



1. The leading order solution 



We first derive the leading order solution for our basic variables. Here, the "leading order" means not the lowest 
order in gradient expansion, but the lowest order of each physical quantity. For example, from Eqs. (|2.10[) , the leading 
order of ip is 0(e°), but the leading one of Ay is 0(e 2 ). 

The O(e ) part of Eqs. dHJ) and are 



V{4>)), 



\K 2 (t) = 3H 2 = 8ttG(- 



K = -3H = 3£P + AttG <y - 2V{4>) 



(CI) 
(C2) 



These equations are indeed the Friedmann equations, and the leading solution is that of a FLRW spacetime. 

Setting j3 l = 0(e 3 ), we substitute the order of magnitude evaluation of the variables shown in Eq. (|2.10|) into the 
Klein-Gordon equation and find 



<t> + 3H <i> + %- = . 

dq> 

d 2 V 

(2)0 + 3H {2) <j) + -j27{2)<t> ~ ((2)« + 3if( 2 )a) - 2(j) {2) a = . 



The Hamiltonian and momentum constraint equations give 

1 



dV 



-0 (2)^ + 0(2)0 +-77 (2)0 



-87rGV 6 M-((2)0)+O(e 5 ): 



where i?y is the Ricci tensor with respect to Yy , and Di is the covariant derivative with respect to 7y . 
The evolution equations for the spatial metric give 



6d t ip 



3H(a-l) = D k /3 k , 



(d t - [3 k d k )% 3 = -2Ay + 'y ik d j p k + % k dtf k - -%d k k + 0(e 4 ) ; 



while the evolution equations for the extrinsic curvature give 



d t A tj + 3HA.J 



2„/,4 [Rid 



a 2 ip 



R] + 0(e 4 ) 



= 3F 2 (2) a + 87rG 



(0 2 + v(4>) 



nil d V 
(2)" + 20(2)0" -77(2) 



First, for ip, we have 



^ = {0) L(x l ) + O(e 2 ), 



(C3) 
(C4) 

(C5) 
(C6) 

(C7) 
(C8) 

(C9) 
(CIO) 

(Cll) 



where /q\L is an arbitrary function of the spatial coordinates. And from Eqs. (|C9|) and (|C11[) together with the 
assumption dtjij — 0(e 2 ), we obtain 



Hi - (O)Jij 



f l3 (x k ) + 0(e 2 ), 



A, 



(2)*ij 



a" 



(0) 



L 4 



a(t')dt'+ {2)C \ ixk) +0(e 4 ) 
n , pi (o) jj *ki 

(2)-«ij + (2)% (0) 



(2)-Ry = f D iDj(fi)L - (0) / ii A( )L 



(0)J 



(0)J 



(2)#fci + (2)-Rfc|) 



—D m LD 0(o) L- 



(0) 



(0) 



L 2 



(C12) 
(C13) 

(C14) 

{f))fi]£>k(p)LD k ( )L , (C15) 
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where t* is an arbitrary reference time, (o)fij is an arbitrary and symmetric tensor of the spatial coordinates, Rij 
is the Ricci tensor with respect to (p)fij> Di is the covariant derivative with respect to ro)fij> A is the Laplacian 
with respect to (o)/y ; and (2)Cy is an arbitrary, symmetric and traceless tensor which depends only on the spatial 
coordinates. 



2. The solution to 0(e 2 ) in gradient expansion 

Now we consider the general solution valid to 0{e 2 ) in gradient expansion. As we have seen in the previous 
subsection, among the basic variables we have introduced, the only quantities whose leading order terms are lower 
than 0(e 2 ) are ip, jij, a, and <ft. Hence what we have to do is to evaluate the next order terms of these variables. 

First let us consider 7y. Substituting Eq. (|C13j) in Eq. (|C8j) . and choosing (3 l — 0(e 3 ), we obtain the general 
solution for 7^ , 

/* df 1 /"* df 1 

^ J a(t")dt" - 2 {2) C ij (x k ) J _ + 0(e 4 ) + {2) f l3 {x k ) , (C16) 

where (2)/ij is an arbitrary and symmetric tensor which depends only on the spatial coordinates. We may absorb it 
into (o)fij without loss of generality. Thus we have 

% = (0)fij(x k ) ~ 2 (2) ^(x fc ) J* f a{t")dt" - 2 (2) C ij (x k ) f + 0(e 4 ) , (C17) 

Here, we note that (o)/y is not completely arbitrary, but its determinant must be unity, det((o)/ij) = 1 + 0(e 4 ). 

To obtain the other variables, we first consider the solution for (2)4 > - We express ( 2 )Q! in terms of (2)0 by using 
Eq. (jUTOl) . We obtain 

W a =j2^^~%^) (C18) 

Inserting this into the field equation (|C4[) . we obtain a closed equation for (2)0) 

, Hj> + 2dV/d0 ■ j>d 2 V/d4> 2 + 2HdV/dcp 

(2)0 ■ (2)0 ■ (2)0 = 0. (C19) 



Although this equation looks difficult to solve at first glance, it turns out that it can be analytically solved. Here, 
we note that because the lowest order scalar field is only a function of time, the above equation is exactly the same 
as the one in the long wavelength limit of linear theory. Then, it is known that there exits a particular solution that 
satisfies 

0(2)0d - ^-(2)0d = 0, (C20) 

where the subscript d is attached since it corresponds to a decaying mode solution in the linear limit. This may be 
integrated easily to give 



(2)0d ocv(t) = exp 



*M (C21) 



a?{t)<}>{t) 



where we have normalized the amplitude so that v = l/a 3 (t„) at t = t*. 

Once we know a particular solution, the other independent solution, u(t), can be found by the use of the Wronskian. 
For two independent solutions u and v of Eq. (|C19|) , the Wronskian, 

W = iiv-vu, (C22) 

satisfies 

W + b(t)W = 0, (C23) 
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where b(t) is the coefficient of (2)4> m Eq. (|C19[) 



H()> + 2 dU/d0 



(C24) 



Thus we obtain 



Then since 



we readily find 



iiv — v u — W (X exp 



b(t)dt 



(a 3 0) 2 



d fu\ uv — vu W 

^2 



dt \v 



W 



ir 



dt' = 



t. « 2 a 3 (*M*) Jt, 



a{t')dt' 



(C25) 



(C26) 



(C27) 



where we have normalized u so that {u/vj = a(t*) at t = t*. Thus, the general solution for ( 2 )0 is given by 

{2) 4> = C(x k )u{t)+D(x k )v{t). (C28) 

Given the general solution for (2)0, the remaining variables <2)Ot and (2)^ are determined as follows. Equation (|C18[) 
gives the solution for ( 2 )a as 



{2)a = ^ 
20 



dU 



20u- — u) C{x k ) + [2(f>v- — v) D{x k ) 



dV 



(j) 3 a 3 



dV 



2acj)+^- I a(t')dt' ) C{x k ) + < ^-D(x k ) 



In terms of this solution, Eq. (|C7|) gives the 0(e 2 ) part of rp as 

^ = ^(0)^(2)", 

where we have set (3 % = 0(e 3 ). Integrating this equation, wc obtain 



ip = (0) L 



1 /"* 

1 + 2 y pja^ 



+ (2)i(^), 



(C29) 



(C30) 



(C31) 



where (2)L(x k ) is an arbitrary spatial function of 0(e 2 ), which we may absorb into t \L without loss of generality. 

Up to now we have not considered the Hamiltonian and momentum constraint equations. The constraint equations 
will relate the quantiti es (gyL , (o)fij, (2)Cij, n)C and (2)D. First let us focus on the Hamiltonian constraint, Eq. (|C5[) . 
Using Eqs. |CTU|) and (fCTlj) . it becomes 



-8 ( o)rAfl j( o)i + (o)f kl (2)Rki(o)L\ = 48^G 



(0) 



L 5 



l dv j. 

{2)9 + -TT(2)<P 



(C32) 



Comparing this equation with Eq. (|C20[) . we see that the right-hand side vanishes for the decaying mode solution v. 
Thus inserting the general solution (|C28p into the above, we obtain 



' [-8 (0) /«A5 j( „)i + { o)f kl (2)Rki io) L] = ~48nG%±C(x k ) . 



a (o)i^ a 
Thus D(x k ) is not constrained by the Hamiltonian, while C(x k ) is determined in terms of ( )£ and (o)/i? as 

-1 



(C33) 



(2) 



c 



487rG (o) L 5 



5 (o) 



f°D t D m L + ( o)/ fei (2)i? fc i(o)iJ + 0(0 



(o) J 

-1 



487rG (o) i 4 
-1 



-(o)/ W [(2)Rk + (2)Rki\ +0(e i ) 



48ttG0* 



i? [L 4 f] + Q(e 4 ) 



(C34) 
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where 0* = <j>(t*) and R [L 4 f~\ is the Ricci scalar of the metric (o)£ 4 (o)/y- 
The 0(e 3 ) part of the momentum constraint (|C6p yields 

& [(o)i 6 ( 2 )^] = -8vrG (0) L 6 ^^ ( (2) C(x fc )) , (C35) 
& [ {0) L\ 2) C tJ ] = -8nG {0) L 6 4>^ ( [2) D(x k )) . (C36) 

The latter equation implies (2)-D is not arbitrary but expressed in terms of iq)L, (o)fij and (2)CV,-, while the former 
equation is found to be consistent with Eq. (|C34[) . This consistency is a result of the Bianchi identities. 
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